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Models of inflation based on axions, which owe their popularity to the robustness against UV 
corrections, have also a very distinct class of signatures. The relevant interactions of the axion are 
a non-perturbative oscillating contribution to the potential and a shift-symmetric coupling to gauge 
fields. We review how these couplings affect the cosmological perturbations via a unified study based 
on the in-in formalism. We then note that, when the inflaton coupling to gauge fields is high enough 
to lead to interesting observational results, the backreaction of the produced gauge quanta on the 
inflaton dynamics becomes relevant during the final stage of inflation, and prolongs its duration 
by about 10 e-foldings. We extend existing results on gravity wave production in these models to 
account for this late inflationary phase. The strong backreaction phase results in an enhancement 
of the gravity wave signal at the interferometer scales. As a consequence, the signal is potentially 
observable at Advanced LIGO/ VIRGO for the most natural duration of inflation in such models. 
Finally, we explicitly compute the axion couplings to gauge fields in string theory construction of 
axion monodromy inflation and identify cases where they can trigger interesting phenomenological 
effects. 
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I. INTRODUCTION 
A. Motivation 

The next few years will bring a large amount of new information about the physics of the very early universe. Cosmic 
Microwave Background (CMB) and Large Scale Structure (LSS) probes will measure the primordial cosmological 
fluctuations with increasing accuracy, and over a widening range of scales. This increasing precision will allow us 
to strongly constrain (or perhaps measure) non-Gaussian statistics of the scalar curvature fluctuations. At the same 
time, gravitational waves (tensor perturbations) will also be probed with substantial improvements in accuracy, by 
CMB experiments and also - at much smaller scales - with gravitational interferometry. 

It is often claimed that the simplest and most natural microscopic models of inflation lead to a rather minimal set 
of observable predictions for future missions. According to the standard lore: 

1. Tensor fluctuations (gravitational waves) are generated entirely by quantum vacuum fluctuations and their spec- 
trum is simply controlled by the Hubble parameter during inflation. One can realistically hope to characterize 
them by just two numbers: the Hubble rate when some pivot scale leaves the horizon, and the tensor spectral 
index, ut- (The spectrum of tensor vacuum fluctuations is detectable only for sufficiently high scale inflation, 
which may be unnatural without invoking symmetries.) This signal is too week to be observed in any of the 
forthcoming gravity waves interferometers pQ. 

2. non-Gaussianity arises entirely due to the (weak) self-interactions of the inflaton field and is undetectably small: 

Jnl ~ 0{e) where e = — is a slow roll parameter and Jnl characterizes the size of the bispectrum f2j-[4]. It 
is often claimed that a detection of primordial non-Gaussianity would be a "smoking gun" signature of some 
non-standard physics in the early universe. 1 

We argue that, even in very natural models of inflation, the situation may be much richer. We present a simple scenario 
where a symmetry renders the high scale of inflation natural, and hence leads to a detectable B mode polarization 
from the tensor vacuum fluctuations. We note that this same symmetry leads also to a series of interesting and, more 
importantly, correlated observables. In particular, we show how detectable resonant and equilateral non-Gaussianities 
may naturally arise, in concert with signals at gravitational wave interferometers (such as Advanced LIGO/VIRGO). 

At the heart of our study is a rather simple and general observation: the inflaton field should couple not only to 
itself, but also to "matter" (non-inflationary) sectors. 2 Indeed, the presence of such couplings would seem to be a 
pre-requisite for successful reheating. Although "matter" couplings are apparently a ubiquitous feature of realistic 



For example: small sound speed [5], higher derivatives [6], special initial conditions [71 [51 ITT| , potentials with sharp features 1 1 21 1 1 31 , 
dissipative effects 1 141 1151 . turning trajectories | 16| . post-inflationary effects (such as preheating 1 171 118| ). etc. 
2 In some cases, "matter" couplings can lead to dangerous loop corrections which spoil the flatness of the inflaton potential, if they are 
too strong. However, this need not be the case in the most robust and appealing inflation models - they rely on some symmetry which 
protects the slow roll parameters and strongly constrains the allowed interactions. 
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models, their relevance for the phenomenology of the primordial cosmological scalar/tensor fluctuations has only 
recently been appreciated. 

In this work, we explore the consequences of "matter" couplings in a particularly natural model: inflation driven by 
a pseudo-scalar axiom In this case, the inflaton enjoys a continuous shift symmetry which, although slightly broken, 
nevertheless protects the flatness of the potential from receiving unacceptably large loop corrections; see [H [TUl [22- 
135] l38l [39] for model building. This construction is especially appealing in the case of high scale inflation - which 
leads to observable tensor modes - since an infinite number of corrections need to be suppressed due to the fact that 
the inflaton travels over a super-Planckian distance [40] . In models of axion inflation, there is generically present a 
coupling to gauge fields of the form 

C int = - — ^F, U F^ , (1) 

where tp is the inflaton, F^ u — d^A^ — d u A^ is the field strength associated with some U(l) gauge field A^, F^ v — 
rf va ^F a p /{2^~g) is the dual field strength (the generalization to non-Abelian groups is straightforward), / is the 
axion decay constant and a is a dimensionless parameter. Notice that, upon integration by parts this is manifestly a 
derivative coupling and therefore does not induce perturbative corrections to the inflaton potential. The coupling 
has important implications for both gravitational waves and also non-Gaussianity. 

It is natural to split up our discussion of the phenomenology of the coupling ([I]) according to the relevant ob- 
servational scales. First, consider CMB/LSS scales. There are two important sources of non-Gaussianity in axion 
inflation: 

1. Couplings of the type ([!]) to non-Abelian gauge groups can induce, via instanton effects, non-perturbative 
oscillatory contributions to the axion potential. 3 These in turn lead to very distinct oscillations in the scalar 
power spectrum and also non-Gaussianity of the resonant type [T3] [213 . 

2. The homogeneous dynamics of the inflaton can lead to nonperturbative production of classical gauge field 
fluctuations due to the coupling ([I]). These produced gauge field fluctuations, in turn, source scalar inflaton 
perturbations by inverse decay: 8 A + 5 A — > 5tp. This effect leads to a large equilateral contribution to the 
bispectrum, as was first realized in |19j . 



In section |III[ we review the computation of both resonance and inverse decay effects using the in-in formalism, 
showing how both effects can coexist in the spectrum and bispectrum. We also generalize previous studies [THJ [20] to 
allow for the possibility of interactions which explicitly break the underlying shift symmetry. (Such interactions are 
actually present in certain string theory constructions; more on this later.) 

Next, let us consider the phenomenology on much smaller scales, generated around 10 to 20 e-foldings before 
the end of inflation. Since the produced gauge field fluctuations contribute to the anisotropic stress tensor, they 
provide an important new source of gravitational waves which is much larger than the usual one generated by vacuum 
fluctuations. Recently Cook and Sorbo [44] have realized that this signal should be visible at interferometers such 
as Advanced LIGO/VIRGO and Einstein Telescope. 4 We continue their analysis by taking into account the effect 
of the strong backreaction regime that necessarily characterizes the late stage of inflation. We show that in the 
regime relevant for observations at interferometers, a copious production of gauge fields leads to about 10 additional 
c-foldings of inflation as compared to the case where the gauge field coupling is absent. As a consequence, the gravity 
wave modes left the horizon when the inflaton field was in a less flat part of the potential (at least, for the simplest 
monomial potentials for which this effect has been studied) and hence a larger population of gauge modes was present. 
This enhances the gravitational waves signal, and makes it detectable for the same range of parameters that lead to 
observable non-Gaussianity, when one assumes the most natural values for the duration of inflation. Specifically, if 
one requires that the large CMB scales left the horizon about 60 e-foldings before the end of inflation - which, as we 
show in subsection II C is the most natural expectation given the coupling ([fj - then backreaction effects enhance 



the gravitational wave signal to a level that can be observed at Advanced LIGO/VIRGO. It is worth pointing out 
that tensor modes produced by inverse decay encode information about the interactions of the inflaton unlike those 
produced by vacuum fluctuations, which depend only on the energy of the background. 

It is interesting to note that the same coupling ([I]) which leads to interesting non-Gaussianity and gravitational 
wave signals will also provide a natural decay channel for the inflaton. We provide a cursory discussion of the physics 



3 The origin of these oscillations is somewhat different in string theory constructions, coming from EDI corrections to the Kahler potential, 
but the end effect is equivalent. 

4 See also |45j for a discussion of the effects of this additional source of gravitational waves on CMB scales. See [46| for a discussion of 
gravitational waves from particle/string production effects in different models. 



4 



of (p) reheating in axion inflation, arguing that the decay of the inflaton is extremely efficient in the observationally 
interesting regime. 

Throughout the first part of this paper, we work in a very general effective field theory context. However, because 
inflationary model-building is sensitive to Ultra- Violet (UV) physics, it makes sense to explore the embedding of QFT 
models of inflation into a more complete framework, such as string theory. In particular we compute the size of the 
coupling ([!]) in a class of axion monodromy inflation in IIB string theory [28] [29l [33] . A model independent coupling is 
present to the gauge fields living on the world volume of the NS5-brane that generates the monodromy. Whether this 
is large enough to trigger the inverse decay phenomenology depends on the size of g s Co (the string coupling times the 
RR zero form). For g s Co < 0(1) as expected in a generic perturbative flux compactification, the coupling ([I]) is too 
small. On the other hand, although it is less common (we found it in .5% of the toy models we considered), nothing 
forbids g s Co > 7 in which case signatures at interferometers and in the CMB would be present. For completeness we 
computed the model dependent coupling to the gauge fields on a D5-branes wrapping the cycle defining the axion. 
This is naturally of the order of the current experimental bound on inverse decay non-Gaussianity, and its precise 
value depends on volume moduli. 

Although we will work in the context of a specific class of models, our key observation - the importance of inflaton- 
matter couplings for non-Gaussian phenomenology - is of a much broader nature. In |47j it was argued that the 
scenario of [19l [20] is couched within a more general class of theories which exhibit the so-called "feeder" mechanism. 
(Another example was explored in detail in [48H5T] .) Indeed, the example interaction ([lj serves to highlight an 
important point: even very simple QFT models of inflation can lead to surprisingly complex dynamics and rich 
phenomenology. 

The paper is organized as follows. In section [IT] we study the homogeneous background dynamics of the model 



(10 1. In section III and IV we study the scalar cosmological perturbations on scales relevant for CMB/LSS and 



interferometers, respectively. Finally, after reviewing in section [V] the string theory construction of axion monodromy 
inflation and the challenges it faces, in section VI we identify cases where the coupling to gauge fields is large enough 
to trigger interesting phenomenological effects. 



B. Effective Field Theory of Axion Inflaton 

To be concrete, let us assume that the effective field theory of the inflaton is characterized by a (mildly broken) 
shift symmetry 

tp —> ip + const . (2) 

Such a symmetry could arise because ip is a Pseudo-Nambu-Goldstone-Boson (PBNG) associated with a global 
symmetry that is spontaneously broken at some scale /. In the context of string theory, the axion arises instead from 
the dimensional reduction of a two-form over an internal two-cycle. In this case, the symmetry ^ has its origin in the 
gauge symmetry of the two-form. Models with an underlying shifty symmetry have attracted a considerable amount 
of interest, partially motivated by the fact that such theories can naturally lead to an observable tensor-to-scalar ratio 

[HH3S1 EH1 EH] - 

The symmetry ([2]), although valid to all orders in perturbation theory, is generically broken by non-perturbative 
effects to a discrete subgroup tp — > tp + 2nf, leading to a periodic contribution to the effective potential of the form 

K P M = A 4 cos(J^ +... (3) 

where / is the axion decay constant, A is a non-perturbatively generated scale and . . . denotes (sub-dominant) higher 
harmonics. The original natural inflation model [2T] exploited the potential ([3| to drive inflation. However, this 
simple scenario is compatible with observation only when / £ 4M p [23], a regime that seems problematic: In the 
case of a PBNG, / > M p would suggest a global symmetry that is broken above the quantum gravity scale, where 
conventional QFT is presumably not valid [231 [55] [S3] ; moreover, values / > M p do not seem possible in a controlled 
limit of string theory [53] . 

In addition to ([3]), one may also incorporate additional ingredients which break the symmetry ^ explicitly at 
tree level and lead to a monodromy: the would-be periodic direction ip is "unwrapped" [28-31 . In many interesting 
scenarios, explicit symmetry breaking ingredients induce a power-law potential of the form 

V sx (<p) = //"V • (4) 



The case p = 1 is typical of string theory models [28l [29] while p = 2 arises from axion/4-form mixing [30] 131 j - 
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Notice that, generically, both contributions ^ and Q may be present: 

V{tp) = V np {tp)+V SI {tp). (5) 

As the suffix sr (slow roll) indicates, it is often assumed that V SI dominates the potential; one then obtains a realization 
of large-field inflation, even for / <C M p . The potential |5| must be constant in the limit that the symmetry |2]) is 
exact. Hence, we expect that the smallness of symmetry breaking effects will protect the slow roll parameters 



from receiving unacceptably large radiative corrections. 

The same symmetry which protects the slow roll parameters also constrains the interactions of tp and therefore 
has implications for non-Gaussianity. In our analysis, a key role will be played by the presence of pseudo-scalar 
interactions ([I]) between the inflaton and the gauge fields. The interaction is consistent with the underlying 
symmetries and, from an effective field theory perspective, it must be included. The strength of the pseudo-scalar 
interaction ([!]) is controlled by the axion decay constant /; absent fine-tuning we do not expect a -C 1. In passing, 
notice that a coupling of the form ([I ) to the visible sector provides a natural decay channel for the inflaton and hence 
may be considered phenomenologica ly desirable. In principle, one can also consider the coupling of the pseudo-scalar 
inflaton to fermions. However, this results in a helicity-suppressed decay rate, which vanishes in the limit of vanishing 
fermion mass. We expect that, if the inflaton decays into significantly lighter fields, the decay into fermions can be 
neglected in comparison to that into gauge fields. 

In the limit that the symmetry ^ is exact, only the coupling ([!]) to gauge fields is permitted (we note that the 
quantity FF is a total derivative so that tp — > tp + c shifts the Lagrangian by a surface term which has no impact on 
the classical equations of motion). However, the explicit symmetry breaking ingredients which lead to the potential 
Q may, in principle, lead to additional couplings. A generic interaction Lagrangian which is compatible with gauge 
symmetry is 

CD _m F 2_m F p. ( 7 ) 

4 4 

The smallness of the symmetry breaking - which is required to ensure slow roll - implies that B{tp) must be nearly 
constant. To quantify the slow variation of the function B{tp), it is useful to define the following parameters 

M v(B'\ 2 „ _ n/r2 B 



£B ^2{b) < (8) 

where the prime denotes derivative with respect to tp. The same logic implies that C(tp) must be close to linear, and 
hence C {tp) should be nearly constant. We therefore introduce an additional set of slow- variation parameters 

(C"\ 2 „ 2 C" 
ec = ^\~a) ' Vc^Ml—. (9) 

(Again here the prime denotes derivative with respect to p.) Throughout this work we assume that ej, |r7j| <C 1 for 
i = V,B,C. 

Putting everything together, and including mild symmetry breaking effects, we are led to consider the following 
effective field theory description of an axion inflaton 

C = - W - fj, i ~ p tp p - A 4 cos ( - ^F 2 - ^?±FF. (10) 



2 v -r, r y f 



The theory ( 10 ) subsumes nearly all known example of axion inflation and generalizes the model of jTHl HO] to 
incorporate also explicit symmetry breaking interactions, parametrized by the coupling functions B{tp), C{tp). 

Throughout this paper, we assume a flat Friedmann-Lemaitre-Robertson- Walker (FLRW) background geometry: 
ds 2 = —dt 2 + a 2 {t)dx? with scale factor a{t). We sometimes use conformal time r = J — . Derivatives with respect to 
cosmic time are denoted as / = dtf and with respect to conformal time as /' = d T f. 



6 



II. BACKGROUND EVOLUTION AND BACKREACTION EFFECTS 



In this Section we study the homogeneous background evolution of the model ( 10 ). We are interested in an inflaton 
potential of the form 



l/M = 14.((^) + A 4 cos(j 



(11) 



We require that V SI is sufficiently flat to drive slow roll inflation, while the oscillatory term may be treated as a small 
modulation. In this Section, we disregard the oscillatory term in (11), which has a s mall impact on the number of 
e- foldings of inflation. 5 Instead, we focus on the effect of the gauge field couplings in (10). We will see shortly that 
the production of gauge fluctuations can lead to important backreaction effects on <fi{t) = (<p(t, x)). These effects 
are the subject of this Section. The discussion is divided in three Subsections. In Subsection |II A| we review the 
production of the gauge quanta due to the motion of the inflaton, generalizing previous results. Some technical details 



of the derivation are relegated to Appendix A. In Subsection II B we show that, if the inflaton-gauge field coupling 
is sufficiently large to lead to observable non-Gaussianity from inverse decay, then the produced quanta affect the 
inflaton dynamics during the last N ~ 20 e-foldings of inflation. This effect produces a change of the predicted scalar 
spectral index n s and tensor-to-scalar ratio r which could be effectively ascribed to a shift of AN ~ 10 of the number 
of e-foldings from horizon crossing for CMB scales to the end of inflation. In Subsection II C we briefly discuss the 
endpoint of inflation, commenting in particular on (p)reheating. 



A. Gauge Field Production 



As first shown in |32j . the homogeneous dynamics of the inflaton leads to an important instability for the gauge 
field, j4 m . To see this effect, we consider the equation of motion for gauge field fluctuations in the background of the 
homogeneous inflaton <f>(t) = (ip(t,x)). We introduce a canonical field variable by the rescaling 6 



Ai(t,x) = VS(0(i))A,(i,x) 



and decompose Ai(t,x) as 



li(T,x) = 



d 3 k 
(2tt) 3 / 2 



e 2 A (k)a A (k)i A (r,k)e lk - x + h.c. 



(12) 



(13) 



where e\ are circular polarization tensors (see |20| for more details) and the annihilation/creation operators of the 
gauge field obey 



a A (k),4(k')l =<5aa«S (3) (k-k') 



The c-number mode functions obey the follow equation of motion 

i±(fc,r) =0, 



92 +k 2 + M 2 A (r)± 2 ^ 



where 



M„ — si: 
2 P B ' 



M 2 a (t) - - 



n (b'^j 



(14) 

(15) 

(16) 
(17) 



5 The oscillatory term in |TTJ is important, however, for the cosmological fluctuations; see Section 

6 In principle this introduces a kinetic mixing with the inflaton, but since we are assuming eg < 1 this is a subleading effect which we 
therefore neglect. 
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and where prime denotes derivative with respect to <f> and we have introduced the slow roll parameter 

ft 

£ *=2lPMl- (18) 



During the conventional inflationary regime we have the usual relations = ey = ~jp (see Eqn. |6|). However, when 
backreaction effects are important these different definitions need not coincide. 

The quantity M^(t), which was not present in [THl [201 132], arises due to the explicit shift symmetry breaking 
interactions Q; this term vanishes when the shift symmetry ([2| is exact and B(ip) = const. Evidently, M\{t) is 
proportional to slow-variation parameters. As we will see shortly, this contribution to the effective frequency of the 



gauge field fluctuations does not have any significant impact on their dynamics. On the other hand, the term in ( 15 ) 
proportional to £ plays a very important role. 

Without loss of generality, we can take £ > 0. Moreover, we will be most interested in £ £ 1. In this case the 
A = + polarization state experiences a tachyonic instability for k/(aH) <, 2£, leading to an exponential growth of 
fluctuations. (The A = — polarization instead remains in its vacuum and can be disregarded.) More specifically, for 
£ = O (1), the gauge field modes are in the vacuum state during most of the sub-horizon regime, and the tachyonic 
instability becomes appreciable only close to horizon crossing. The growth ceases when the mode becomes much 
greater than the horizon, and then the energy density in that mode decreases due to the expansion of the universe. 
Therefore, at any given moment during inflation, only gauge modes of comparable size to the horizon are relevant. 

Typically, the parameter £ oc 4>/H increases during inflation. However, the time variation of £ appears only at 
higher order in e, rj. We can therefore treat £ as an adiabatically evolving parameter: we treat it as a constant 
during the interesting part of the evolution of each mode (since |A£| <§; £ over At — few x H^ 1 ). However, when 
comparing two well separated stages of inflation, we take into account that £ can be different between them, leading 
to a different amount of production. We note that £ instead varies significantly after inflation, so that we do not 
expect that treating £ as a constant provides a good approximation to the gauge field production during this stage. 



In Appendix A we derive an analytic expression for the properly- normalized solutions of (15 1, in the approximation 
that e^, £b and £ may be treated as constant (slow roll approximation). We show that the solutions admit the simple 
representation: 



We note the exponential enhancement, e 71 ^, in (19 1, which reflects the underlying tachyonic instability. Equation (19) 
coincides (up to a change in normalization) with the result that was first derived in |32] for the case B — 1, C oc <p. 
As first discussed in [19], this approximation is valid in the region (8£) _1 ;$ - fcr < 2^ that accounts for most of 
the power in the produced fluctuations (we stress that this phase space is nontrivial for £ > 0(1), which we assume 
throughout; the production of gauge fluctuations is uninterestingly small for £ < 1). The result ( fl9| ) highlights our 
previous claim that the explicit symmetry breaking interactions Q do not significantly impact the production of 
gauge fluctuations in axion inflation. See Appendix A for more details. 



The modes (19) are real- valued (up to an irrelevant constant phase); this leads to an important simplification when 
we compute non-Gaussian correlators. The fact that the c-number modes are real-valued implies that the produced 
gauge field fluctuations are commuting variables, to a very good approximation. That is, we have 

[d r MT,x),Aj(Ty)] «o, (20) 

where it is understood that only the produced fluctuations with —kr < 2£ are relevant (the modes with — kr 3> 2£ 



remain in their vacuum and their effect is renormalized away). Equation (20) reflects the essentially classical nature 
of the tachyonic production of gauge fluctuations (we note the analogy with the decoherence and semiclassicality of 
the standard cosmological perturbations due to the expansion of the universe [53]). 



B. Classical Evolution 



In this subsection, we study the backreaction of the produced gauge field fluctuations ( 19 1 on the homogeneous 
dynamics of the inflaton <j){t). The smallness of the parameters ^ and ^ implies that -B(co]should be very close to 
a constant (that can be absorbed into the normalization of the gauge field) and C(ip) should be very close to linear. 
Hence, for simplicity, we focus on the case 



C(<p) 



a <j) 
2fH 



(21) 
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for the remainder of this section. We also focus the discussion to the string theory axion monodromy inflationary 
potential (see Section VI) 



(22) 



During inflation we have ip S> <p c so that V(ip) « /k 3 |</j|. (The parameter ip c is, however, relevant during reheating.) 

We introduce the physical "electric" and "magnetic" fields as (despite the terminology, the gauge field need not 
correspond to the Standard Model photon) 



B 



1 



;Vx4, 



E 



1 



A'. 



(23) 



The gauge field fluctuations ( 19 1 are produced at the expense of the kinetic energy of the homogeneous inflaton 



4>(t), introducing a new source of dissipation into the equation of motion of 4>. Moreover, these produced fluctuations 
contribute to the energy density of the universe and thus modify Friedmann equation. These effects are encoded in 
the mean field relations 



3H<j> + V' = j(E-B), 



3H 2 = 



Ml 



I • 



F + V + -(E 2 + B 2 



(24) 



and, using (19), one finds 



Therefore QHCUl] 



(E-B)~ -2.4 • 1CT 4 |J e 27r « 



,E 2 + B 2 



TTi 

1.4-10- 4 ^e 



2?r£ 



(25) 



£ 3/2 c ?r£ ^ — _^ ne giigibi e backreaction on cf> eq. 
H 2 

^-3/2 e 7rf ^ _^ negligible backreaction on Friedmann eq. 

H 



(26) 



During the standard slow roll regime 



,FM P « HM p 



and the first condition is the more stringent than the 



second. If we take the standard result for the spectrum of primordial perturbations, P. 



1/2 



c 



7 (2W> 



5 • 10" 



this condition reads £ < 4.7. 

In pi|l [20] it was shown that £ ~ 2.5 leads to the standard result for the power spectrum (up to a subdominant 
correction due to the modes generated by inverse decay), and to observable (but not yet ruled out) non-Gaussianity 
from the inverse decay of the gauge quanta into inflaton quanta. We see that backreaction is indeed negligible for 
such value. However, we note that £ typically increases during inflation. For concreteness consider a linear inflaton 
potential (such is the case for monodromy), then £ oc l/<j) during the standard slow roll regime. Let us assume that 
£ = 2.5 at |0| ~ HMp, which is the value of the inflaton about 60 e- foldings before the end of inflation, under the 
standard slow roll assumption. We then find that £ = 4.7 is reached at \4>\ ~ 5.9 M p . This happens before inflation 
ends. As discussed in [32] . the production of gauge quanta results in an additional friction on the inflaton motion 
that, in the current context, prolongs the duration of inflation. 

We evolved numerically the two equations ( 24 ) , keeping into account the backreaction term in the equation for 
the scalar field. We still disregarded the energy density of the gauge quanta in the Friedmann equation; in this 
way this equation remains a simple algebraic equation for H ; the right panel of Figure [2] shows that this is a good 
approximation all the way to the end of inflation. We choose parameters such that, initially, we are in the standard 
slow roll regime adopted in [19j[20j. Therefore, the result for the power spectrum, and for /nl obtained in [19) 120]. 
apply (see Subsection III B ) . We note that the coupling between the inflaton and the gauge field can be given in terms 



of the initial values of £ and 



<t>CMl 



2.5 * 7 =5 



">CMB | 



(27) 



Here and in the following, the suffix "CMB" denotes the value of quantities when the large scale CMB perturbations 
left the horizon (when this statement needs to be quantified, we refer to the WMAP pivot scale [56 k = 0.002 Mpc -1 ). 
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FIG. 1: Evolution of the inflaton as a function of the number of e-foldings to the end of inflation, starting from \4>cmb\ = 
— 9.9 M p , with (red solid line) and without (green dashed line) the coupling to gauge fields. For the first line, the strength of 
the inflaton-gauge field coupling is chosen so to lead to observable non-Gaussianity from inverse decay. For the second line, we 
have shifted the number of e-foldings to make manifest that the two evolutions coincide at early times. 

We obtained \4>cmb\ — 9.9 M p ; this is smaller than the value (\4>cmb\ — 10.9 M p ) leading to 60 e-foldings of 
inflation without gauge production, confirming that the backreaction of the produced quanta increases the amount of 
inflation. 

In figure [l] we show the evolution of the inflaton field as a function of the number of e-foldings to the end of inflation 
for £cmb = 2.5 (red solid curve) and for £ = (green dashed line), i.e the standard slow-roll case. The backreaction of 
the produced quanta on the background evolution becomes noticible during the last ~ 25 e-foldings of inflation, while 
it is negligible at earlier times. The two trajectories reach <f> = at different times, showing that the backreaction 
increases the duration of inflation by about 10 e-foldings. 




FIG. 2: Left panel: Friction terms in the equation of motion for (f>. Right panel: relative strength of the energy density of the 
produced quanta; this term is neglected in the numerical evolution of the background equations. 

This change in behavior during the last ~ 25 e-foldings of inflation is also visible in the left panel of figure [2j where 
we show the evolution of the two friction terms in the inflaton equation as a function of the number of e-foldings 
to the end of inflation. The standard Hubble friction controls the earliest stages, but the system gradually evolves 
towards a regime in which the backreaction of the produced gauge quanta dominates the evolution. Namely, the 
system approaches the strong backreaction regime studied in |32) . Let us stress that in the our case the observable 
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cosmological fluctuations were produced during the phase of standard slow roll inflation (and therefore, the results 
of [19j [20] apply), while the strong backreaction regime of [32] here is reached only in the last O(10) e-foldings of 
inflation. In the right panel of figure [2] we show the relative contribution to the energy density of the gauge quanta to 
the total energy density. This term was neglected in the numerical evolution of the system (24). The Figure confirms 
that this is a valid assumption. 



C. Reheating 

We conclude this Section with a brief comment on reheating in our model. It is interesting that in this model the 
same coupling (pFF which lead to non-Gaussianity will also provide a natural decay channel for the inflaton. We 
expect reheating into the gauge bosons to be very efficient whenever particle production effects are important. 
This is illustrated in the right panel of Fig. [2j There we see that, already by the end of inflation, the energy density 
in gauge field fluctuations is comparable to that of the inflaton. At this point it may be expected that the condensate 
4>{t) is efficiently destroyed and thermal equilibrium may set in rather quickly. 



A detailed account of the dynamics of (p)reheating in the model ( 10 ) goes beyond the scope of this paper. Here we 
briefly discuss the perturbative decay of an axion inflaton, anticipating that such an analysis might be viewed as a 
lower bound on the efficiency of reheating. The perturbative decay rate into gauge fields associated with the coupling 
is 

a m% 

T ^ AA = 64^ ■ (28) 



The inflaton mass is obtained fr om the potential V SI . For string theory monodromy inflation models this is given by ( 22 ) 
where, as we discuss in Section VI ip c = O (10 _1 ) M p (we note that, indeed, V^ r ~ M 3 |0l during inflation). This leads 



to ~ 5 . 10~ 5 M p . Taking a/f as in (27), and equating T v ^aa = H, we obtain the estimate T r h = O (lO 12 ) GeV 
on the reheating temperature of the thermal bath generated by the perturbative decay. As mentioned above, for the 
values of the coupling we are interested in, the energy density of the gauge quanta becomes comparable to that of the 
inflaton already at the end of inflation. We expect that this will lead to additional backreaction effects beyond the 



mean field backreaction of (E ■ B) and (E 2 + B 2 ), which is the only effect that enters in (24). Specifically, we expect 
strong rescattering effects that should lead to a much quicker thermalization than the perturbative rate suggests. As 
for preheating of scalars, we believe that it would be worthwhile to study these processes through lattice simulations. 

It is worth mentioning that we have implicitly assumed that the inflaton-gauge interaction in (JlJ is the strongest 
interaction of the inflaton to any gauge field. Pseudoscalars have also a natural coupling to fermions AC = 
_Cm^p (j)^^^ where C is a model dependent quantity which, in the spirit of an effective field theory, should not 
be expected to be <C 1. The associated decay into fermions is helicity suppressed, ~ (C '/ '/) m^m 2 / (8ir) for 

m,,p <C 771.0. This is parametrically smaller than the decay rate into gauge fields if -C m^. 

The detailed history of reheating is necessary to relate the size of any given mode to the number of c-foldings N 
before the end of inflation at which this mode exited the horizon. This relation is given by |55j : 

where the suffices CMB, "end" , "rh" , and denote horizon crossing during inflation, the end of inflation, the conclusion 
of reheating, and the present epoch, respectively. For £cmb — 2.5, and \4>cmb\ = 9.9M p , this relation evaluates to 

AT(fc)~61.7-log T-~> 1 / ■ (30) 

V ; B 0.002 Mpc- 1 3 p V 4 V ; 

(Notice that, due to the logarithmic dependence, this relation will change only very slightly if £cmb and </>cmb differ 
from the specific values considered here.) We have normalized the wavenumber k to the WMAP pivot scale [5T>] . 

The dependence on p r h in this relation is due to the assumption that, between the end of inflation and the completion 
of reheating, the massive inflaton field is coherently oscillating about the minimum of its potential, resulting in a matter 
dominated stage. We see that, in general, a longer period of matte r do mination during inflation results in a lower 



value of N. Above, we have seen that the perturbative decay rate (28) leads to T r h ~ 10 GeV. Inserting this in 



(30), we obtain TV ~ 59 for the CMB pivot scale. This however assumes matter domination of the inflaton oscillations 
before the perturbative decay takes place. We however see that, already by the end of inflation, the gauge fields 
carry a non- negligible fraction of the total energy density. One may therefore expect that the equation of state of the 
universe will have an intermediate value between that of radiation and matter [57] . This would lead to a number of 
c-foldings between 59 and 62 for the CMB pivot scale. 
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III. FLUCTUATIONS AT CMB/LSS SCALES 



In this Section, we study the scalar cosmological fluctuations in the model (10) on CMB/LSS scales, corresponding 
roughly to scales which left the horizon ~ 55 — 60 e-foldings before the end of inflation. In section [TTJ we have 
seen that backreaction effects are negligible during this regime. Even though backreaction effects are negligible, the 
production of gauge field fluctuations still has an important impact on the observable cosmological fluctuations, via 
inverse decay processes [THEQ]. Here, we re- visit the computation of [19), [20] using the in-in formalism and allowing 
for the possibility of explicit symmetry breaking interactions. We also account for the possibility of resonance effects, 
showing how these "add up" with the fluctuations from inverse decay. 

This Section is divided into four Subsections. In Subsection MI Al we derive the interaction Hamiltonian for the 
theory (10), accounting for all relevant effects. Technical details are provided in Appendix B. This effective action is 



the key input for a computation of the n-point correlation functions using the in-in formalism. In Subsection |III B| 
we briefly outline the computation of the spectrum and bispectrum using the in-in method. Because the technical 



details are quite involved, we expand on this analysis in Appendices C and D. Finally, in Subsection IIIC we study 
observational consequences and delineate interesting regions of parameter space. 



Effective Action for the Perturbations and In-in Calculation 



We are interested in a potential of the form ([11]), where the oscillatory term can be treated as a small correction 

(31) 



and backreaction effects are negligible. It is useful to introduce a parameter 

A 4 



/TO*)' 

where cj)^ is 4>(t) evaluated at the moment t+ when some relevant pivot scale fc* left the horizon. If we assume b < 1 
so that ( |TT| ) is monotonic, then the data require [55] b <C 1 so that a perturbative expansion in b is justified. We will 
assume this regime throughout the paper. The solution for the background evolution of the inflaton, <fr(t) = (<p(t,x)), 
can be computed to leading order in b, with the result [41j 

, M , M 3b/ 2 . (Mt)\ ^ 
cj){t) = 0o(t) sin — — , (32) 



where = —H/H 2 is a slow roll parameter, 7 evaluated at t*. We have also introduced <fio(t), which is the solution for 
the background scalar field in the absence of modulations, that is for 6 = 0. (Obviously, 4>o(t) depends on the specific 
choice of V ST .) The slow roll parameters also admit an expansion in b <C 1. In particular, we can write e = —H/H 2 
as: 

e = e + ei H , (33) 

ei = -36^-V2^cos (^j . (34) 

We are interested in the n-point correlation functions of the variable 



C(t,x) = -^M*,x)- (35) 



These may be computed using the in-in formula 

(C* Ck 2 • • • Ck„ (t)> = ]T H) N dh dt 2 • ■ • / dt 



tpt—i 



JV=0 " u " u "° 

x([[[C kl Ck 2 ---Ck>),i^i)],#/fe)]--- ,Hj(t N )}), (36) 



In this Section we work in the regime where backreaction effects are negligible and it is not important to distinguish between ej,, ey 
and —H/H 2 . 
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where Hj(t) = — J d 3 x£i is the interaction Hamiltonian and, with some abuse of notation, it is understood that the 
fluctuation modes on the right hand side correspond to the free theory solutions. 

We choose the free theory to correspond to the quadratic action for £, neglecting resonance effects (setting 6 = 0). 
To leading order in slow roll parameters, the free theory Lagrangian is given by 



Co = M'a 6 e 



-(VC) S 



The solutions of ( 37 ) are well-known. The q- number Fourier mode is written as 

Ck(r)=6(k)u fc (T) + ftt(-kK(r), 
where the annihilation/creation operators are 

[&(k),&t( k ')] =6^ (k-k') , 

and the c-number modes are given by 

H 2 e~ lkT 



(1 + ikr) 



(37) 



(38) 



(39) 



(40) 



up to an irrelevant overall phase. 

Now we turn our attention to the leading interaction terms. These were derived in Appendix B. The result is 



Hi{t) 



* CfT^FaP + My £l (c 2 ^(VC) 2 ) + ^(2e ) 3 / 2 A 4 ^ sin f M 



6 



P 



(41) 



This result is derived neglecting inhomogeneities of the metric, which has been show to be a good approximation 
both in the case where inverse decay effects dominate [5D], and also when resonance effects dominate [35] . The first 



term in (41 1 leads to inverse decay effects. The contribution to the 2-point and 3-point correlation functions from 



this term is studied in Appendix C. There we show explicitly that the in- in formula (36) gives the same result as the 
method employed in [THIHO]- The last two terms in (41 1, on the other hand, lead to resonance effects . The relevant 
contributions to the 2-point and 3-point functions are computed in Appendix D. Here we notice these two distinct 
physical effects "add up" in a very simple way, at the level of the n-point correlation functions. 



B. The Spectrum and Bispectrum 



The 2-point correlation function defines the power spectrum 



2?r 



{&Cv(T)) = —P c {k)5^ (k + k') 



(42) 



The N = term in (36) gives the usual (nearly) scale-invariant result. There are two important corrections to this. 
At N = 1 there is a contribution due to the quadratic interaction terms in (41 ); these gives an oscillatory modulation 



of the power spectrum that was accounted for in [29, 41j using a different approach. See Appendix D for a derivation. 
At N = 2 there is an important contribution from the pseudo-scalar coupling which influences the normalization 
of the spectrum; this was accounted for in |19l I20j and is re-derived in Appendix C. Accounting for all important 
contributions, we can write the power spectrum in the form 



P c( fc) = A 2 (- 



This is the same form that was analyzed in |29l 141 j . however, the parameters now have the meaning 

A 2 = V[l+Vf 2 (0e 4 ^] , 
2tt/ 



5rii = 3&* 



1/2 



(43) 

(44) 
(45) 
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Recall that P 1 / 2 = H 2 /(2n0). The function /2(C) was computed in [2TJ] and is given by 

MO : 



3-10" 5 ^- 5 - 4 if 2 < £ < 3; 
7.5 ■ 10~ 5 ^ 6 if £ > 1. 



(this function is plotted in Figure 1 of [20]). 
The total bispectrum is defined by 



(C kl Ck 2 Ck 3 ) =B(h) (ki + k 2 + k 3 



(46) 



(47) 



This quantity receives two important contributions. At N = 1 we have a contribution due to the cubic self-interaction 
term in (41 ); this gives a resonant contribution that was studied in [29] |4l] 02] . See Appendix D for more details. At 
N = 2 we have unimportant corrections to this result. The leading contribution from the pseudo-scalar interaction 
in (41 1 arises at N = 3 and gives an equilateral contribution to the bispectrum that was first derived in |19) : see 



Appendix C. Summing up all important contributions, leads to 



^inv.dcc \k% ) H~ -^rcs(^i) 

9(2tt) 5 / 2 
10 



P 2 dk)fTL AcC {0T^{h) 



-(2tt) 5 / 2 P c 2 (A;) 



res 
NL 



h 2 h 2 h 2 



<\>K 
f 



k~ 



where the nonlinearity parameters are 



I NL 



p 2 (k) 



rrcs 
JNL 



The equilateral shape template is given by 

1 1 



3&*V2tt ( J2e*M< 



f 



3/2 



2 1 

+ 



1,31,3 Up pp h 2 h 2 h 2 k,h 2 lc 3 
12 1 3 2 3 1 2 3 I 2 3 



(5 — perms) 



and the function /s(£) is 



/ 3 (0 = 



7.4- lO- 8 ^ 81 if 2 < ^ < 3; 
2.8 -10- 7 ^ 9 if^> l. 



(48) 

(49) 
(50) 

(51) 
(52) 



(this function is plotted in Figure 2 of [213]). 

To be precise, the shape of the non-gaussianity from inverse decay does not coincide exactly with the equilateral 
template (see Figure 6 of [20]). The cosine (defined in [58]) between the precise shape and the equilateral template 
evaluates to 0.94 [2U], so that the equilateral shape can be used as a good estimate of the actual result. A positive 
detection of non-gaussianity based on the equilateral template will justify a more precise analysis where the exact 
inverse decay shape is taken into account. 

The results of this subsection illustrate two important points. First, we have seen explicitly that the in-in method 
reproduces the same answer as the formalism developed in [TH] HO] for inverse decay effects. This serves as a useful 
consistency check and helps to verify the equivalence of two of the most popular methods for computing cosmological 
correlation functions. Second, we have shown that inverse decay effects and resonance effects "add up" in a simple 
way, both at the level of the spectrum and at the level of the bispectrum. 

It is interesting to note that the two different effects which we have discussed - resonance and inverse decay - lead 
to very different patterns of higher order correlation functions 47 . This difference may be relevant for observables 
that are sensitive to the pattern of moments, such as the statistics of rare objects or Minkowski functionals. 



C. Observational Consequences at large scales 



Models of axion inflation have a rich phenomenology. The presence of a periodic modulation in the potential (11) 
gives rise to resonant effects, while the inflaton coupling to gauge fields results in inflaton perturbations from the 
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inverse decay of the produced quanta. The resonant signatures were studied in ref. [13, 29, 4"lTfl5] while the inverse 
decay signatures were studied in [TU [20]. For definiteness, we discuss the resonant effects only for the linear case here, 
V SI = fi 3 (p during inflation. (This case can be derived from string theory [28, 29 a and will be discussed in more detail in 
sections [V] and [VI}) The parameter /i is fixed by matching the observed power spectrum normalization. There remain 
3 key parameters which determine the cosmological observables: the axion scale /, the amplitude A of the periodic 
term in the inflaton potential (equivalently, we can use (31) to replace A with b), and the coupling of the inflaton to 
gauge fields a. At a given choice of /, the magnitude of resonance effects are controlled by b and the magnitude of 
inverse decay effects are controlled by a. 



© -2 




loglO(f) 



FIG. 3: Parameter space for axion monodromy inflation with V B r — fJ- 3 <P- The "TT(res)" regions give the (one- and two-sigma) 
likelihood contours for the temperature 2-point function. The "TTT(res)" regions give / rcs = 200,20,2. The region to the left 
of each horizontal line is ruled out by observational constraints on .f^' 1 , for values of the coupling a = 0.04, 0.2, 1 (see eq. 
The Figure has been adapted from |29j by adding the exclusion regions from . 



The phenomenological constrains on these parameters are summarized in Figure [3] The "TT(res)" regions show 
one- and two-sigma likelihood contours for the power spectrum, with the darker regions at large b being disfavoured 
by the too large oscillations in the power spectrum induce by resonant effect. The "TTT(res)" regions give the value 
of resonant nongaussianity with (from top left to bottom right) / ros > 200, 20, 2, respectively. The exact bounds on 
the size of resonant non-Gaussianity from current and future experiments are not yet known, but research in this 
direction is in progress. The horizontal lines indicate the value of equilateral nongaussianity due to inverse decay. For 
a given choice of a, the region to the left of the horizontal red line is ruled out by the 95 % C.L. WMAP [55] bound 
f e *™ < 266. Recall that the generic expectation from effective field theory is a — 0(1). We see that, in this case, 
observational constraints on equilateral non-Gaussianity can rule out much of the parameter space where resonant 
non-Gaussianity is large, but still allow for detectable oscillations in the power spectrum. To have a reasonably large 
resonant bispectrum without a large inverse decay signal requires a much smaller than 0(1). For such values of the 
pseudo-scalar coupling, we may have a richness of possible signature in the CMB with oscillations in the spectrum 
and bispectrum, along with equilateral non-Gaussianity. 

In Figures [4] we present instead results both for linear and quadratic V SI during inflation. In the left panel we show 
the value of inverse decay non gaussianity as a function of //a. We recall that the inflaton-gauge field coupling is 
inversely proportional to f/a, see eq. ([I]), which explains why the resonant decay non-gaussianity is a decreasing 
function of this quantity. The rapid decrease observed in the Figure is due to the fact that the effect is exponentially 
sensitive to the coupling [THl HO] • In the right panel we show instead the phenomenological predictions for the spectral 
tilt n s and the tensor to scalar ratio r as the inflaton-axion coupling ranges in the values shown in Figure [4] The 
predictions are compared with the 68 % C.L. and 95 % C.L. limits given in [56 . 
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FIG. 4: The left panel shows the inverse decay non-gaussianity as a function of the inverse of the inflaton-gauge field coupling 
(see eq. {!])). The value shown corresponds to modes that left the horizon 60 e- foldings before the end of inflation. The WMAP 
limit is obtained from the 95 % C.L. WMAP [56] bound on equilateral nongaussianity, and from the fact that the shape of 
inverse decay non-gaussianity has a cosine of 0.94 with the equilateral template The right panel shows values of n s and r 
for the same range of inflaton-gauge field couplings. These values are compared with the 68 % C.L. and 95 % C.L. limits given 
in [5(3]. Both panels show results for linear (p — 1) and quadratic (p — 2) slow-roll infiaton potentials. 

The results shown in Figure [4] improve over the existing literature [20] by taking into account the backreaction effects 
of the produced quanta on the background evolution. As we have discussed in Subsection |IIB[ for inflaton-gauge field 
coupling leading to interesting non-gaussianity, backreaction is negligible when the scales relevant for CMB (and LSS) 
observations left the horizon. However, it becomes relevant at the last stages of inflation. The main backreaction 
effect is the slow down of the infiaton due to its loss of kinetic energy into gauge quanta. This in turns increases the 
duration of inflation with respect to the case of no production. As a consequence, the number of e- foldings Ncmb 
inside inflation at which any CMB scale left the horizon is a function of both the initial value of the infiaton 4>cmb 
and of the coupling of the infiaton to gauge fields a/ f. If we fix the value of Ncmb associated to a given scale (this 
request can be done once we know the reheating history after inflation), the increased of the duration of inflation 
needs to be compensated by a value of 4>cmb which would lead to a smaller value of Ncmb in the absence of the 
coupling. For the monomial slow roll potentials considered here, this means a smaller value of <f>cMB- Therefore, 
Figure [4] has been obtained for a progressively smaller value of 4>cmb as the coupling a/ f increases. 

The bottom-right end of each of the two lines shown in the right panel of Figure [4] corresponds to the standard 
chaotic inflation results obtained in absence of coupling. Such values are obtained as long as roughly £cmb 1 
(corresponding to f/a^, 0.046M p in the linear case, and f/aZ, 0.065A/ p in the quadratic case). The upper-left 
ends of the two curves correspond instead to the maximal couplings compatible with the CMB non-gaussianity limit 
(namely, we exclude £cmb ~ 2.66, corresponding to f/a < 0.019M p in the linear case, and f/a < 0.027M p in the 
quadratic case). As the coupling increases towards the maximal allowed value, the corresponding decrease in 4>cmb 
implies a less flat potential when the CMB modes leave the horizon. This in turns results in increased values for 
1 — n s and r (this is the same reason behind the lines in the n s — r plot as a function of the number of e-foldings 
in the standard case, cf. Figure 20 of [SB]). In the present case, there is an additional effect on r at large couplings. 
Namely, the gauge quanta give a greater contribution to the scalar power spectrum than to the tensor power spectrum 
[T9], [20] . Specifically, the standard relation r ~ 16 e is modified into r ~ 7e 2 in the limit in which the inverse decay 
perturbations dominate over the vacuum one [20]. As shown in [19l [20], the inverse decay non-gaussianity limit 
forces one to be in the regime in which the power is dominated by the vacuum mode. For the maximal possible 
couplings, the infiaton perturbations from gauge modes contribute to about 10% of the scalar power spectrum, while 
the corresponding contribution to the tensor power is about 3 • 10~ 4 for linear infiaton potential, and about 6 • 10~ 4 
for quadratic infiaton potential. This is the cause of the decrease of r at the left end of the curves shown in the right 
panel of Figure [4] 
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FIG. 5: Q.gw h 2 as function of the frequency /, for N = 60 e-foldings of observable inflation, a linear slow roll inflaton potential, 
and £omb ~ 0, 2.33, 2.66 (the value of £ when the large scale CMB modes left the horizon). For reference we also show the 
expected sensitivity of LISA, Advanced LIGO/VIRGO and Einstein Telescope (at their most sensitive frequency). 




FIG. 6: Region in the {Ncmb,£cmb} plane (values assumed by these quantities when the large scale CMB modes left the 
horizon) for which the gravity wave signal is detectable at Advanced LIGO/VIRGO and Einstein Telescope. The left and right 
panel refer to a linear and quadratic inflaton potential, respectively. 



IV. GRAVITATIONAL WAVES AT INTERFEROMETERS 



In Section III we discussed the observable cosmological fluctuations on CMB/LSS scales. Such scales left the horizon 
roughly 55 to 60 e-foldings before the end of inflation, during the phase where backreaction effects are negligible. In 
this section, we instead study scalar and tensor fluctuations on much smaller scales. These modes left the horizon 
closer to the end of inflation, when backreaction effects start to play an important role in determining the evolution 
of the homogeneous background, </>(t) and H(t). Our main results are summarized in figure [6j where we show that 
Advanced LIGO /VIRGO could detect a stochastic background of gravitational waves from inflation for £cmb as small 
as 2.33 (equivalent to f/(M p a) < 0.021) in the case of a linear inflaton potential, and as small as 2.23 (equivalent to 
f/(M p a) < 0.031) in the case of a quadratic potential. 
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A. The Gravitational Wave Signal 



The inverse decay of tachyonic gauge fields does not only affect scalar but also tensor perturbations. The power 
spectrum has been computed in [19], and further studied in |2Ql|45]. Defining gij — a 2 (8ij + hij) with hij transverse 
and traceless, ha = dihij = 0, and its Fourier transform as 



hij (»^3 ^) 



d 3 k 
(2tt) 3 / 2 



^ H-ij, r {k) [h r (k,t)a^ + h.c] 



(53) 



r=L/R 



one finds the tensor power spectra 



P, 



L/R 



h 3 



H 2 



IT 2 M 2 



m2 f 

-j^p-Jh,L/R(!;)> 



4tt£ 



(54) 



for the left and right polarization respectively. Here the functions fh,L/R are defined by 



fh,L/I 



d 3 q* 



g* | A; - q* 



(1 ± cos 6) 2 ^1 — cos 8 ± ^/l — 2^ cos 6 + q\ 



(27T) 

dxyfx [sin x — x cos x] 



16(1-2(7, cos6 + q 2 ) 

-2 % /2|^[ v /^;+ v / |fc-9.ll 



(55) 



where and 9 are, respectively, the magnitude of the (dimensionless) integration momentum q*, and the angle 
between this vector and the momentum k of the mode. For details, and for the numerical evaluation of these 
functions, see [213] , 

In order to compare our results with the sensitivity of interferometers, it is convenient to introduce the fractional 
energy density per logarithmic wavenumber interval 



1 



3H 2 M 2 aiogfc ' 



where the label refers to quantities evaluated nowadays and 

Ml 



Pgw,o 



(KjKj + d k hijd k hij) 



(56) 



(57) 



To evaluate pcw,o we need to evolve the primordial signal ( 54 ) after horizon re-entering until today. Using the tensor 
transfer function (see e.g. |36U37j ) for modes that re-entered the horizon during radiation domination, as it is the case 
for the modes relevant for interferometers, one finds 



n 



w 



R,0 



24 



(Pl + Pr 



(58) 



where f^o = Pr,o/^HqM 2 ~ 8.6 • 10~ 5 refers to radiation nowadays (one needs to include also the neutrinos as if 
they were still relativistic today). It is useful to plot Qgw,o as a function of frequency / = k/2n, which is related to 
the number of efoldings by 



N — Ni 



CMB 



In- 



= In 



0.002Mpc 

kcMB 



In 



0.002Mpc _1 
/ 



0.002Mpc 



-r - 44.9 - In , 
- 1 10 2 Hz 



(59) 



The resulting signal as a function of frequency is presented in Figure [5] for the case of a linear inflaton potential - 
which is an excellent approximation to the potential for axion monodromy |28j - and for three reference values of the 
inflaton-gauge field coupling. Specifically we choose £cmb — 2.66, corresponding to roughly the highest value allowed 
by the CMB limits on equilateral non-gaussianity, £cmb = 2.33, corresponding to roughly the lowest possible value 
that can be detected at Advanced LIGO/VIRGO, and £,cmb — 0, corresponding to standard vacuum fluctuations. 
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The results presented in the left panel of Figure [5] have been obtained assuming Nqmb — 60. A lower duration 
of inflation leads to a greater signal, as we show in the right panel 8 . The fact that the signal grows with £cmb is 



due to the exponential term in (54), while its growth with decreasing N C mb is slightly more subtle. For fixed £cmb, 
decreasing Nqmb affects the scales relevant for interferometers (N = Nqmb ~ 44.9) in two ways: it decreases H, 
since these scales are now closer to the end of inflation and it increases £, since the 4>/H is larger. The latter effect 
ends up increasing the total gravitational signal due to the exponential enhancement. 

Finally, let us comment on how the consistent account of the strong backreaction phase affects the results of 
[44] . where this phase was neglected. One effect can be seen in Figure [5] for a fixed £cmb 0, where one can 
distinguish three phases. From left to right one sees the vacuum contribution, then the fast growth of the inverse 
decay contribution and finally the slow down of the growth when the backreaction becomes sizable. A second effect, 
which is clearly visible by comparing the two curves for different £,cmb is that for larger £cmb the whole signal is 
effectively shifted toward lower frequencies (to the left) due to the larger number of additional efoldings of strong 
backreaction at the end of inflation. The two effects act in opposite direction, but since the latter is stronger than the 
former, the net result of the strong backreaction phase is an increase of the signal. This increase is crucial in order 
to make the natural window Nqmb ^59 — 62 potentially visible at Advanced LIGO/VIRGO. 

B. Perturbations in the Regime of Strong Backreaction 

So far, we have discussed the tensor perturbations on small scales, relevant for interferometers. One should also 
study scalar cosmological fluctuations at such scales. Although these are not directly observable (besides the possible 
generation of primordial black holes), it is nevertheless important to verify that scalar fluctuations are perturbatively 
small, as a consistency check. For scales which left the horizon ^10 — 20 e-foldings before the end of inflation, the 
calculation presented in Section |III| is not valid, since backreaction effects no longer have a negligible impact on the 
homogeneous dynamics. Let us provide a brief, heuristic analysis of the scalar fluctuations in the regime of strong 
backreaction, in order to justify the validity of perturbation theory (which was implicit in our previous analysis of 
gravitational waves). According to Ref. [32], the inflaton perturbations obey an equation of the form 

V 2 

5ti + 3/3H5cp ^5ip + m 2 5cp = J , (60) 

a z 

where the effective coefficient of friction, /3, and the source term, J, are given by 

0=1-27^^-^, J=-\e-B-{E-B)}. (61) 
J SH<p J L J 



We have assumed that £, 4> > and (25) shows that (E ■ B) < 0, so the effect of particle production is always to 



increase the effective friction in (60), as would be expected on physical grounds. Equations similar to (60) have 
also been obtained in other scenarios where particle production effects contribute an important source of friction; see 
[T4"l [IB] for example. 

When backreaction effects on the homogeneous dynamics are negligible we have j\(E ■ B}\ <C 3i?|</>|; see Section 
|nj In this case /3 = 1 and ( |60| coincides with the equation studied in \nj\ [20]. Instead, when backreaction effects 
are important /? 3> 1 is possible. We can give a heuristic estimate of the size of scalar fluctuations in the regime of 
strong backreaction as follows. Since most of the interesting production occurs near horizon crossing, we estimate 



the derivatives as d ~ H. In this case, and for /3 ^> 1, (60) gives: 3(3H 2 (Sip 2 ) 1 / 2 ~ (J 2 ) 1 / 2 . We further estimate 



(J 2 ) 1 / 2 ~ j\(E ■ B)\ and £ ~ — ¥-8<p. The variance of fluctuations is roughly 



(C 2 )~O(10- 2 )l, (62) 



which is consistent with the analysis of [32 . Equation ( 62 ) suggests that curvature fluctuations on small scales are 



much larger than on CMB scales, although a perturbative analysis is still justified. (A useful benchmark value to 
keep in mind is £ ~ 5 — 6 near the end of inflation for the case where non-Gaussianity is observable, but not yet ruled 
out.) 



As we discussed in subsection II C for £ = O (1), having Nqmb < 60 requires however a non-minimal reheating history. 
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V. AXION MONODROMY INFLATION 



We start this section by reviewing a construction of axion monodromy inflation in string theory following [281 129j 
and a generalization proposed in [33j . We then conclude by summarizing and commenting on the main open issues 
in the current implementation of this class of models. 



A. Review of the Model 



Consider an 03/07 orientifold Calabi-Yau compactification of Type IIB string theory. The relevant four-dimensional 

upon dimensional reduction. Following [29] . we choose the normalization 



uij = (2tt) V<5/ (63) 

for a basis £/ of the homology H2(Y,Z) and luj of the cohomology H 2 (Y,Z), where Y is the Calabi-Yau space over 
which we compactify. With an ansatz for the ten-dimensional RR two-form 

C 2 = ±-c{x)u , (64) 

Z7T 

for some base two-cycle U), we get a four-dimensional axion c(x) with periodicity 10 2ir. The axion decay constant for 



CIS 11 



f = 9s 
Ml 12V £ (2^) 2 



(27r) 6 a' 3 



S^-P^' (65) 



where Ve is the six-volume of Y measured with the Einstein frame metric. In order to break the axion shift-symmetry 



(the details of the resulting potential will be reviewed in section VI I, we consider a spacetime filling NS5-brane 
wrapping the cycle X that defines the axion c. If this cycle is far away from any O-plane, the local physics is given 
by the oriented Type IIB string theory, since the orientifold projection relates the fields at some point x to the fields 
at some other far away image point x' and so does not lead to any local constraint. Tadpole cancellation requires 
an additional anti-NS5-brane wrapping a homologous two-cycle. A metastable configuration can be obtained if we 
assume that £ belongs to a one parameter family of cycles that stretches between two warped throats and brane and 
anti-brane are located in different throat. Since the tension of the branes is reduced by warping, there is a potential 
barrier for the branes to annihilate. Inflation takes place as the C2 flux over S slowly decreases to minimize the 
NS5-brane (an anti-brane) tension. A series of consistency checks and bounds have been discussed in [2"51 129] . 

To conclude, we want to mention a slight modification of the above construction described in |33j that alleviate 
some of the difficulties encounter in this model. Imagine that for some reason it is hard to compute or control the 
monodromy potential generated by the introduction of NS5-branes (some concrete difficulties will be discussed in the 
next subsection). Inflation can still take place if there is another axion (in general one expects h^i axions, i.e. as many 
axions as orientifold-odd two-cycles) which interacts with the one discussed above via a non-perturbative effect (which 
are not forbidden by any apparent symmetry). The resulting two field potential (plotted in polar coordinates) closely 
resemble the Inferno as described by Dante, which was taken as the name of the model. If there is a mild hierarchy 
between the two axion decay constants (the precise amount depends on the details of the monodromy potential) a 
long stage of slow-roll inflation can be obtained. 



9 -B-type axions, coming from the NSNS two-form, receive a large mass from the non-perturbative corrections used in the moduli stabi- 
lization a la KKLT. Although one might very well conceive a scenario in which moduli are stabilized by perturbative corrections, this 
is harder to construct in detail. 

This can be seen e.g. via S-duality starting from the world-sheet coupling J B / (2irct'). Note that this choice differs from that in |28l . 
where the axion periodicity was (2-k) 2 . 

1 Let us check that this is compatible with e.g. (3.5) of |59l and (2.15) of 28 . The volume scaling L -4 agrees sinc e v a ~ L 2 E . For the 
scaling with g 3 one has to transform from the string-frame volume of 1281 1601 to the Einstein-frame volume of |65| i.e. L 4 = g s L E , so 
that (f/Mp) 2 ~ g 2 /Lj, = g s /L E . The fact that 1281 1601 use a different 4D Einstein frame from the one we use here doesn't matter since 
/ and M p scale the same under a constant red-shift of the metric so that the ratio is invariant. 
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B. Issues with Concrete Realizations 

Before proceeding to the study of (f>FF in the string realization of axion monodromy, we would like to stress that 
there are still several challenges that the explicit constructions face. Although addressing these challenges is beyond 
the scope of the present paper, we list some of them below. 

Explicit solution. One should provide an explicit solution with the required geometry (proposals were give in |28j ) 
and solve the ten-dimensional equations of motion in the presence of warping, fluxes and an NS5/anti-NS5-brane pair. 
Among other things, it would then be interesting to compute the exact warped axion decay constant and the masses 
of NS5-brane modes. 

Backreaction on four-cycles. The C2 flux (i.e. the inflaton) on £ induces and effective (anti-)D3-brane charge on the 
(anti-)NS5-brane, which changes the warp fact, which in turns modifies all warped volumes. Since four-cycle volumes 
are stabilized, this chain of interactions leads to an effective potential for the inflaton. In the single axion monodromy 
model this effect can be suppressed if the NS5/anti-NS5-brane pair behaves effectively as a dipole, i.e. when the 
distance between them is much smaller than the distance to the closest stack of D7-branes. In Dante's Inferno this 
effect is smaller to begin with [33] for two reasons. First, the flux on the two-cycle wrapped by the NS5-brane is 
suppressed by the ratio of the two axion decay constants, which is assumed to be small, and second the inflaton 
direction is almost perpendicular to the direction for which the backreaction induces a potential. 

Light KK-modes. As pointed out in [55], flux along the compact direction of brane decreases the masses of the 
Kaluza-Klein modes. This can be understood by considering the T-dual picture where flux becomes tilting of a lower 
dimensional brane, which is therefore longer that in the fluxless case. In the single axion-monodromy case, these 
lighter masses are very close to the Hubble scale, while they can be made parametrically heavier in the Dante's 
Inferno scenario. 

NS5/D7-brane intersection. 12 Just on topological grounds, any two-cycle should intersect its dual four-cycle at 
least at one point. In the Type IIB compactifications where we understand how to stabilize the moduli (based on 
|61| ) some non-perturbative corrections arc needed to stabilize the Kahlcr moduli (four- and two-cycles). In the single 
axion monodromy model [28j this should arise via gaugino-condensation on a stack of spacetime filling D7-branes 
wrapped on a four-cycle for each homology class. The effects of the resulting NS5/D7-brane intersection on both 
the moduli stabilization and on the inflaton potential should be studied. Perturbative Kahler stabilization could in 
principle solve this issue but it is computationally much harder to make explicit. 

In the rest of the paper we will compute certain local quantities that should be rather independent on the way the 
above issues are addressed. 



VI. COUPLINGS FROM STRING THEORY 



We devote this section to the explicit computation of the bosonic inflaton couplings induced by the presence of an 
NS5 and a D5-brane and their phenomenological consequences. 



A. Model Independent Couplings: NS5-brane 

The action for the NS5-brane can be obtained via S-duality from that of a D5-brane as discussed in Appendix E. 
In the following we focus on the potential for c and its couplings to the gauge field living on the world- volume of the 
NS5-brane. Performing an S-duality transformation on the action for D5-branes one finds 

s%M = -«V 5 / ^ 4 v /= 5^v /e0|r|2 ^ 2 + (27rc)2 ( x + (2 T^p^) ' (66) 



?C5 _ 3 M5 (2W) d /• E , A ^ _ 1 



- SlHT^i FAFcC Q = 2( ^ )W / 1 F,-f., 
where a subscript E indicates Einstein frame and r = Co + ie~^ . 



This issues was raised in private communications among R. Flauger, L. McAllister, T. Wrase and one of the authors 
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Defining F pa = F^e^ 1 " 7 / (2-^/— <7_e), one can rewrite this action using the notation of section IB Including also the 
axion kinetic term, we have 



S = I d 4 x^/-g E 
where the coupling functions are 

C(<p) 



4 r r ^ 



C(cp) 



1 



{dp) 2 - V{y) 



B(<p) = 



(2^) 2 / 
(a'f 



4"2 V9E,29s\T\ 



(68) 

(69) 
(70) 



After imposing that the potential has a vanishingly small cosmological constant, the inflaton potential takes the form 

(71) 

where 



M„ 



(2tt) 2 V e 



1+G05, 

P (2vr) 2 ^2 P 



(72) 



In the last equality we assume an isotropic compactification for a rough estimate. During inflation ip 3> p c and the 
potential can be approximated as V((p) ~ /i 3 M- (Notice that, for such large field values, we also have B(ip) cx V(ip).) 



B. NS5-brane Phenomenology 



As discussed in section II A the presence of a tachyonic gauge mode depends on the parameter 

C.,h 



2 M ^ 



(73) 



To estimate this we assume that 60 e-foldings before the end of inflation <fi ~ 10M p 3> /, hence we can neglect the 
first term in the square root in the definition of B. Our estimate is therefore 



M 2 



(74) 



Notice that all powers of / cancel. Now the two possible regimes are Co <C l/g s or the opposite Cq l/g s . In each 
case one finds 



M 



2 ^C g s ~ C o5s 10- 2 « 10- 2 

Ml o So 

^ (C o5s ) 3 ~ (C o5s ) 3 10- 2 , 



(75) 
(76) 



respectively. To have nontrivial production of gauge field fluctuations requires £ >, 0(1). Evidently, this is possible 
only when Co > 7/g s , while no relevant effect is expected from the coupling to the NS5-brane degrees of freedom for 
Co < 7/g s . Larger values of Co seem less natural but not impossible. Allowing for fine tuning, Co can be made large, 
as one can see in the explicit examples of [52] where eq. (5.7) gives the value of Co at the minimum. As a figure of 
merit, we found that sampling random values for the fluxes considered in [52] about .5% of the cases had Co > 70. 
The phenomenology of the regime Co > 7/g s is completely analogous to the one we reviewed in section III The inverse 
decay proceeds through the CFF coupling, while BF 2 gives only sublcading corrections. Potential observables are: 
non-Gaussianity in the CMB, running of the spectral index and gravitational waves at interferometers. 
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C. Model Dependent Couplings: D5-brane 

One could assume that there is also a spacetime filling D5-brane wrapping a two-cycle 13 which contains the two- 
cycle £ that defines the axion c. This brane is not essential for the construction, but for completeness we include a 
study of its effects. In order for the D5-brane not to complicate the construction outlined in section [V] and not to 
heavily affect the inflationary dynamics, it must be stabilized at a different location from the NS5-brane. This could 
arise because it wraps a distinct but overlapping cycle, or because of a particular configuration of fluxes. 

In this case, a coupling of c to the gauge theory living on the world-volume of the D5-brane is induced by the 
Chern-Simon part of the action. One finds 



S D tts(W) d - / FAFc, 
where of course this F is different from the one living on the NS5-brane. The kinetic term for F reads 



4(27r)3 ff; 



/ d 4 xy/-g E F E , 



(77) 



(78) 



where vf is the volume (measured in the string frame) in units of y a! of the two-cycle wrapped by the D5-brane. 
Canonically normalizing one obtains 



S 



cfixyf-gE 



:Fl 



Oi(j> 

4/ 



FeFe 



(79) 



with 



14 



2ng s Zn^/gs 



i'2 



(80) 



where now v 2 is in Einstein frame, and / given by |65| 



D. D5-brane Phenomenology 

Again the phenomenology is dictated by the size of £. Focusing on CMB scales we find 

C - # = K\f^ - 2.8 (81) 

For the toy model of an isotropic compactification the powers of the volume cancel exactly leaving £ remarkably close 
to but slightly above the bound imposed by non-Gaussianity. This means that if a D5-brane is present in the same 
cycle as the axion, the phenomenology discussed in this paper is a natural outcome. A very explicit construction is 
needed to evaluate the volume factors which could easily move £ above or below its experimental bound. 

Another way of discussing the result ( |80| ) is by looking at figure [3j A natural size of a is of order .1 to .01 which 
makes the inverse decay effect of comparable order as the one induced by oscillations (whose overall size of course 
depends on the independent parameter b). 



VII. CONCLUSIONS 



In this paper, we have studied a very simple and natural class of inflation models wherein ip is a pseudo-scalar 
axion. We have found that this construction leads to a surprisingly rich array of correlated observables for both scalar 



Again wo assume that in case an orientifold projection has been performed, the two-cycle in question is far away from O-planes so the 
that the local physics is described by Type IIB string theory. 

Let us compare this with (6.7) of |60| . They define L 2 as the string- frame volume measured in units of \Z%~ko! as can be seen from their 
(6.5)-(6.6). This accounts for the 1-k difference. The factor of g a on the other hand is mysterious for me. They should have a gj 1 in 
the kinetic term in (6.5), since the DBI action comes from a disk amplitude (and they are in string frame). In fact they do have the 
right factor in (3.14), but not in (6.5). So we disagree here. . . 
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and tensor cosmological perturbations on CMB scales and also at much smaller scales. For natural values of the axion 
decay constant - roughly / ~ 10~ 2 M p - these signals may be detectable in the not-too-distant future. 

On CMB/LSS scales, we have argued that a detectable non-Gaussian signature arises naturally. Non-Gaussianity 
in axion inflation can feature equilateral and resonant type bispectra or a superposition of the two. We have reviewed 
these effects using the in-in formalism and generalized previous results to allow for more general classes of gauge field 
couplings (which explicitly break the underlying shift symmetry). Since most axion models are of the large field type, 
we also expect a large tensor-to-scalar ratio. 

We have found that axion inflation may lead also to detectable gravitational wave signatures at Advanced 
LIGO/VIRGO scales as realized in [44 . We extended that analysis by including the backreaction effects on the 
late time inflaton dynamics that necessarily arise for couplings leading to observational effects. In particular, we 
have noticed that such effects can increase the number of e-foldings by ~ 10, which has important consequences for 
observables. Specifically, it increases the values of 1 — n s and r at CMB scales, and of Slcwh 2 at interferometer scales. 
The gravity wave signal is correlated with the large-scale non-Gaussianity and, if both effects were detected, it would 
provide a remarkable confirmation of this scenario. 

We have also noticed that the same pseudo-scalar coupling which leads to equilateral non-Gaussianity and small- 
scale gravitational waves will also provide a natural decay channel for the inflaton. Indeed, we have seen that reheating 
is likely to be extremely efficient in the observationally interesting region of parameter space. 

Most of our phenomenology relies on the pseudo-scalar coupling <pFF, which is constrained only by naturalness at 
the effective field theory level. In order to address UV sensitivity we have computed this coupling in specific string 
theory constructions of axion monodromy. We have found that such couplings must be present to gauge fields living 
on the NS5 brane. However, the coupling is large enough to lead to interesting phenomenology only when the moduli 
stabilization is such that g s Co > 7, which can be achieved roughly at the cost of a .5% fine tuning. For completeness 
we noticed that, remarkably, the coupling to a putative D5-brane gauge fields would naturally take the value currently 
constrained by observations. 

Before concluding we would like to mention a few interesting directions for future research. It would be interesting 
to: 

• Extend to smaller scales the study of the phenomenological limits due to inverse decay. The running of the 
spectrum is larger than for the vacuum contribution and could receive bounds competitive with those from CMB 
non-Gaussianity from the inclusion of high I data of experiments like ACT iflj or SPT [72 . Stronger bounds 
could also emerge by combining non-gaussianity at both CMB and LSS scales, taking into account the specific 
evolution of of £ (and of /nl) betweeen the two scales. 

• Find the observational bound on resonant non-Gaussianity. 

• Study observational constraints on scales in between CMB and interferometers. Although probes of primordial 
perturbations at these scales might be less sensitive, and more model dependent (for instance, interesting limits 
can be obtained if one assumes that the dark matter consists on neutralinos; see |73j for a summary of limits at 
various scales) the signal is predicted to grow monotonically as one approaches smaller scales (still remaining 
perturbative) . 

• Forecast the constraints that various interferometers could put on the parameters of the model, in case of a 
detection. 

• Provide a detailed account of (p)reheating in string theory axion monodromy, along the lines of |63H67j . 

• Study in detail the mechanism that keeps perturbation theory under control and determine how large the 
tensor-to-scalar ratio can become in the regime of strong backreaction. 

We hope to address some of these questions in future research. 
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Appendix A: Production of Gauge Fluctuations 



In this appendix we study the tachyonic production of gauge field fluctuations in the theory (10). This analysis 
generalizes the previous work [19, 20, 32 to account for generic coupling functions B(ip), C(tp). We assume that the 
slow variation parameters (JsJ) and ^ are small. 

Working in Coulomb gauge, the equation of motion of the gauge field in the background of the slowly-rolling inflaton 
condensate takes the form 



A" 



Bdy 1 VA B dtp 



(V x A)i = . 



We introduce the re-scaled field 



ii(t,x) = y/B(<f>(t))Ai(t,x), 



and decompose A4 into circular polarization modes, as in equation (13). The equation of motion is: 

d 2 



dr 2 



k z + M\{t)± 



2k£ 



A ± (k,r)=0, 



where 



\J e 4> e B . 

o sig: 



C" 



n [B 1 ^ , 



where = 2 h 2 m 2 ' wri ich may differ from ey (see Eqn. 



2 M p — sign 



6) or —H/H 2 in the regime of strong backreaction. 



(82) 
(83) 
(84) 

(85) 
(86) 



In equation (85) we work to leading order in slow- variation parameters. In the following, for brevity, we assume 



that 4>C' > and B'4> > 0. One can show that 



(87) 



This shows that, when the slow-variation constraints are respected, we can treat £ as a constant. In this case, equation 



(84) has the form of Whittaker's equation 

r d 2 1 A 1/4- fj, 2 ^ 



dz 2 4 



W Xifl (z) =0. 



The gauge mode solutions may therefore be expressed in terms of Whittaker's function. The correctly normalized 
solution is 



A+(K t) = e^ 2 W- K ^+2ikr) , 



(89) 
(90) 



where B+ is B(<f>(t)) evaluated at horizon crossing for cosmologically interesting modes. Equation (89) can be put in 
a more illuminating form by using the asymptotic expansion 



W. 



. w(;) ,(^" < t -,v* 



(91) 



which is valid for large |k| and when Iul(k) > 0. Equation (91 1 leads immediately to the following limiting behaviour 
for pi): 
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Formally, this expression is only valid for £ —> oo. However, one may verify numerically that it provide a good fit 
even for £ ~ 2 — 3, at least for modes with (S^)^ 1 ;$ — fcr <, 2£ which account for most of the power in produced 
fluctuations. 

Equation (92 1 differs slightly from the result that was obtained in [TDJUD] for the case B = 1 and C oc ip. The key 



"new" feature encoded in this equation is the fact that, if B -C 1 during inflation, then the value of £ is effectively 
increased. In this case, one can find a much greater production of gauge quanta, since the amplitude of A; is 
exponentially sensitive to £. 



Notice that equation (92 1 will fail to adequately describe the dynamics on sufficiently large scales. This expression 
results from neglecting the last term in square braces of (881, an assumption that breaks down for — fcr < 0(e/£). 



From (82 1, we can see that the "new" interaction associated with B(tp) actually dominates the dynamics of in the 
limit — kr — > 0. Hence, we should verify that this interaction does not lead to any interesting particle production on 
super-horizon scales, that would not be accounted for by the asymptotic expansion (92). Neglecting gradients, the 
equation of motion ( 82 ) becomes 



A, 



1 dB 
H+ B^t 



Ai^Q. 



Suppressing the vector index, the solution is 



A = A* + A* dt 



,a* B* 
a B 



(93) 



(94) 



Our assumption about slow-variation of B(t) means that (aB) 1 should be a decreasing function of t and the integral 
will rapidly converge. Assuming that B varies slowly enough to pull it out of the integration we have 



B H \ a 



const . 



(95) 



This shows explicitly that the coupling B(<p) does not lead to any interesting particle production on super-horizon 
scales. 



Appendix B: Effective Action for the Perturbations 



In this Appendix, we derive the leading interactions terms in the effective action for the curvature perturbation. 
Following 19, [2D] and [42j we neglect metric perturbations. The action for 5<p, including terms up to third order in 
fluctuations, is given by 



dtd x 



dtd 3 xa 3 



A 4 



cos 



(M 5 



6/3 



■ sin 



(96) 



Here we have introduced the notation m 2 = V^. and we write the zeroth order homogeneous solution 4>o(t) in the 
trigonometric functions, since the difference <j>{t) — (f>o(t) is subleading in the small parameter b = jyr', see 
Moreover, we have disregarded a contribution proportional to V^" , which is usually subleading in a slow roll expansion 
and is actually vanishing for string theory axion monodromy inflation. 



Our goal is to re- write the action ( 96 1 using 



C(*,x) 



H 



Sip(t, x) 



Sip 



2eM p ( 



(97) 



and identify the leading interaction terms. Let us focus, first, on the quadratic terms on the first line of (96). Using 
( 97 1 is is straightforward to derive the result 



S D Mz / dtd 4 xa A e 



i 



(VC) 2 - £l 2 C 



(98) 



where we have defined 



n 2 = - 2 



A 4 



m~ — ttt cos 
/ 



-H z 



H 



3e + 36 + 2e 2 + 3e5 + S 2 



(99) 



2G 



and introduced the parameter 5 = H/(2HH). In the case without resonance effects (when 6 = 0) the last term in 
(98) would rep rese nt a tiny and unimportant correction. However, it is important to notice that it is </>, and not 
(f>o, appears in (971, hence there is rapid oscillatory time dependence implicit in the parameters e, S appearing in il 2 
which could, in principle, lead to interesting effects. To extract the oscillatory terms from f2 2 , we expand the slow 
roll parameters as 

e = e + ei H , S = 5 Q + Si H (100) 

The subscript i indicates that a quantity is computed at i-th order in b. The leading corrections are given by [5TJ 22] 



Si 



ni fc\ f (Mt) 



-3b sin 



p 

Mt) 
f 



(101) 
(102) 



cancels at leading order. Hence, we are left with 



which shows that S 3> e. Using ( 101 1 and ( 102 1 it is straightforward to check that oscillatory contribution to f2 2 

(103) 



S 2 w M 2 I dtiPxafe 



j3„,„3. 



1 



(vcr 



which gives the same linearized equation of motion that was studied in 4 

C" + 2aiJ(l + <5)C'-V 2 C«0. 



(104) 



We take the "free theory" to correspond to ( 
contributions as perturbative "interactions" . T 



rat is, we identify the free theory Lagrangian as 



C Q = M 2 a 3 e 

The 0(b) correction to the quadratic Lagrangian gives 



103), neglecting resonance effects (setting 6 = 0), and treat all other 

(105) 



1 



C'--a(VC)' 



Hi D Af 2 a 3 ei 



c 2 



:(VC)^ 



(106) 



There is also a cubic interaction arising from the last term on the first line of ( 96 1 . To leading order in b this gives 

,3 



/:^y(2 £ o) 3/2 A 4 ^si 



f 



(107) 



Next, we turn our attention to the gauge field interactions on the second line of (96). At leading order we have 

Ci D a 3 B((j) ) 



(108) 



Notice that, when particle production effects are relevant, we always have ( ^, 1 > ^/eeg. Thus, the first term in 
( 108 1 always dominates in the observationally interesting regime. We can further simplify the result by making a field 
redefinition An — > B^ 1 / 2 A u and neglecting derivatives of B : which is justified to leading order in the slow-variation 
parameters (fSl) - Thus, we arrive at 



where e^" 13 is the Minkowski-space Levi-Civita symbol, with rj 0123 = +1. 

Adding up all the relevant terms, we arrive at the following result interaction Lagrangian 



(109) 



Ci = M 2 o 3 ei 



Q 2 



1 



(vc) 5 



+ ^(2 £ o) 3/2 A 4 ^sin 



C 3 + y^F^F a0 



(110) 



Equation (110) is the main result of this Appendix. 
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Appendix C: Inverse Decay Effects with the In-In Formalism 



In this Appendix we study inverse decay effects using the in-in formalism, showing that this reproduces exactly the 
results which were derived in 19, 20 using a different method. The relevant term in the interaction Lagrangian (110) 
is 



(Ill) 



This coincides with the interaction that was studied in [19, 20 (only the meaning of the constant £ is changed due to the 
presence of non-trivial coupling functions B(tp) and C((p)). The interaction Hamiltonian is given by Hj = — j d 3 xC] 
and may be written in the form 



H} nv 



dec / 



(t) = V^oM p / AC-q(^qW , 



to leading order. Here the "source" term is defined by 



2 v / 2^ M p J (2tt) 3 / 



d 3 3 



-ik-x 



F^F, 



(F^F 



(112) 



(113) 



Notice that (Jq(t)) = by construction. 

The quantity J q has the same meaning as in j!9l 120) : it is the source term in the equation of motion for Sip that 
corresponds to the production of inflaton fluctuations by inverse decay. Defining the canonical field variable 



Qk(t) 



the equation of motion becomes 



dl + k 2 



a 
a 



Qk(t) = J k (t) + 



(114) 



(115) 



where • • • denotes additional corrections, arising from resonance effects and self-interactions, which will not concern us 
here. We are interested in the particular solution of ( 115 ), which corresponds to the inflaton fluctuations generated by 
inverse decay and may be highly non-Gaussian. Working at the level of the equation of motion, it is straightforward 
to see that the inverse decay contribution to the n-point correlation functions of the curvature perturbation may be 
written as 



(CkiCks • • • COIj 



inv.dcc 



i 



n 



dr. 



Gkj (t, Tj) 

a(r) 



X (Jk 1 (Tl)Jk 2 (7- 2 ) ' ' ' Jk„(T„)) 



(116) 



where Gfe(r, t') is the retarded Green function for equation (115); see 20; for details. Here we show that the same 
result may be derived also from the in-in formalism. 

Focusing on the interaction term (112) the in-in formula (36) is 

ft fti ftpf — i 



(c kl Ck 2 ---Ck»> = jt(-i) N fdh r dt 2 --- f 

N=0 J J J 



(It 



N 



([[[Ck 1 Ck a -"Ck.(T),flr- dec (<i)] ,^r- dec (t 2 )]--- ,Hr- dec (t N )]) 



CO 

E 

N=0 



(-iV2eoM p ) 



N 



N 

n 



u-i 



dti / d 3 qi 



x ([[[Ck^k.-'-Ck^r),^^)^^)],^^^)^^^)]--- ,C-q N (t2v)J qN (tiv)]), (H7) 

where to = t in the product over integrals. A crucial simplification arises from noting that the "source" terms J qi {Ti) 
may be treated as commuting variables, to a very good approximation. This follows from the discussion in Subsection 
II A| see equation (20) in particular. This allows us to pull the factors of J qi (ri) out of the nested commutator in 



(117). Thus, we arrive at a simplified formula 

oo 

(Ck X Ck 2 -Ck») = ^Hv^Mp 



N 



N=0 

x ([[[CkxCka 



n 



dU / d qiJqiin) 



■ Ck„ (r), C- qi (*i)] , C- q2 (*a)] • • ■ , C- qN (t N )}} 



(118) 
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To evaluate the nested commutators on the last line of (1181 we can use the formula 



[CkiCrO.Ckafa)] = 



1 1 



[Qkx{n),Qk a {n)] 



-j G kl (T 1 ,T 2 ) 

2e M2 a(Ti)a(r 2 ) 



6& ( kl + k a ) , 



(119) 



where the second equality is valid only for t\ > r%. 

Let us first consider the 2-point function; equation (118) with n = 2. The N = term is the summation just gives 
the usual (nearly) scale-invariant power spectrum from the quantum vacuum fluctuations. Inverse decay effects, on 
the other hand, are encoded in the N > 1 terms. Using (119) it is straightforward to verify that only the N = 2 
contribution to the nested commutators gives a non-vanishing contribution. Explicit evaluation gives 

i \ 2 

(CkiCk 2 W} linv.de 



-1 



r) / dT 1 dT 2 G kl (T,T 1 )G k2 (T,T 2 ) (JkiCn)^^)) 
p J J — oo 



(120) 



in precise agreement with (116). 

Next, we turn our attention to the 3-point function. This time the N = term is trivial, since the free theory 
modes are gaussian. It may be readily verified that only the N = 3 term in the summation (118) survives; all other 
commutators are zero. Hence, we find 



(CkiCksCkaCT))!; 



.dec 



-l 

T^Mp 



3 „T 



dT 1 dT 2 d,T 3 G kl (r, Ti)G k2 (r, T 2 )G ka (t, t 3 ) ( J kl (n) J k2 (t 2 ) Jkafo)) (121) 



which, again, agrees with (116). The analysis of this Appendix proves that the formalism developed in [TH1 US] is 
equivalent to the in-in method. 



Appendix D: Resonance Effects with the In-In Formalism 

In this Appendix we study resonance effects [13l [4TJ [42] on the 2-point and 3-point function, using the in-in 
formalism. These effects are generated by the first two terms in the interaction Lagrangian (110): 



C 2 -^(VC) 2 



(122) 



Let us first consider the 2-point function. The relevant term in the interaction Hamiltonian is 



D - / d 6 xa 6 e 1 



1 
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ft z 



(123) 



The N = 1 term in the in-in formula ( 36 ) gives 



(Cki Ck 2 (t)> DiAf„ 2 / dt^xe^t) 



(124) 



Recall that the oscillatory time dependence of ei is given by (101). It is straightforward (but non-trivial) to evaluate 
(fl24|. We find 
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(125) 



where x = —k\T. We can write the homogeneous solution as <fio = cj) k + y/2eoM p \nx; see [H]. Following [3T], we 
introduce the quantity 



= -3bi- 



f 



dx'e 2lx ' cos 
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f f 



lux' 



(126) 
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In this case the total 2-point function, accounting also for the N = term in (36), can be written as 



2ir 2 



(C kl Ck 2 (r)> = -jg-PcMSW (k 1+ k 2 ) 



where 



C 



(k) = v( 



= V 1 + 2Re 



,(-) 



(127) 



(128) 



This coincides with the result of [41] that was obtained using a different method. 

Next, we turn our attention to the 3-point function. The relevant term in the interaction Hamiltonian is 



Hj(t)D- / d 3 ^" r (2 eo ) 3/2 A 4 -f si 



6 " P 



Mt) 

f 



c 5 



(129) 



This interaction was taken into account using the in-in method in [42] , where it was shown that the bispectrum agrees 
with [25] . 



Appendix E: S-duality 



S-duality is part of the SL(2,R) symmetry of the classical action Type IIB (broken to SL(2,Z) at the quantum 
level). The 10D Einstein frame metric Ge is related to the string frame metric G s by Ge,mn = £~^^ 2 G s .mn- It is 
Ge that is invariant under SX(2,R) and not G s . In the 10D Einstein frame, SX(2,R) acts as [55] 



ar + b 



H' 3 



d c 
b a 



H 3 
F, 



(130) 



with ad — be — 1, where r = Co + ie & and all the other fields being invariant. S-duality corresponds then to the 
choice a = d = and b = —c = ±1, leading to the transformation r' = — 1/r or 



G °- CI 



C 



Cl + e-^ 



The D-brane action in Einstein frame is 6i 

Sdbi — —^5 
Then by S-duality one finds the NS5-brane action [TP] 

Sdbi - -Ms J d & ^\r\J -d,A 
= -Ms / d 6 C 



y rf 6 Ce _0 y / -det(e^/ 2 Gf d + B ind + 27ra'F) , 
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(131) 
(132) 

(133) 
(134) 

(135) 



where F = dAi is the field strength of the vector field on the worldvolume of the NS5-brane. If we assume Co = 
and compactify to four-dimensions we find 15 



S D -n 5 a' J d 4 ^- s /-g E ,i\flf 



(136) 



Note that here we have made a different choice of the four-dimensional Einstein frame with respect to [28] ■ There the four-dimensional 
metric is the string metric and the extra factors of g B are reabsorbed in the definition of the Planck constant (which is hence also different 
from what we have here). The advantage of the present choice is that the SL(2,Z) transformations are simpler in the ten-dimensional 
Einstein frame. On the other hand, one has to remember to go to Einstein frame when using the DBI action! 
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where l\a! is the volume of the two-cycle measured with the Einstein metric. One can check that the tension of the 
NS5-brane (in string frame) has an extra factor of g~ x with respect to the tension of a D5-brane as expected since 
the two 5-branes are related by S-duality. 
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